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会 摘要 ”太阳 高 能 粒子 (SEP) 事 件 是 一 类 重要 的 空间 天 气 灾害 性 事件 ,如 能 准确 预报 SEP. 事件 ,人 们 便 可 以 采取 必 

要 的 防护 措施 ,保障 卫星 . 星 载 设备 以 及 航天 员 的 安全 , 尽 可 能 地 降低 经 济 损失 . 因此 ,其 数值 预报 研究 在 空间 天 气 

预报 研究 中 占有 很 重要 的 地 位 . SEP 事件 中 的 高 能 粒子 在 不 同 的 时 间 尺 度 内 被 炊 斑 过 程 或 者 CME 驱动 的 激 波 加 

F= 速 ,并 且 在 被 扰动 后 的 行星 际 太 阳 风 中 传输 ,这 些 过 程 都 紧 紧 依赖 于 太阳 风 背 景 场 .因此 获取 更 加 接近 物理 真实 的 

太阳 风 背 景 场 是 模拟 SEP 事件 的 重要 部 分 ,也 是 提高 SEP 物理 模式 的 关键 因素 之 一 .我 们 目前 的 工作 基于 张 明 等 

PERI SEP 在 行星 际 空间 传播 的 模型 ,尝试 将 Parker 太阳 风速 度 解 及 WIND 飞船 观测 的 磁场 实时 数据 融和 人 模型 

中 ,研究 不 同 的 太阳 风速 度 以 及 真实 磁场 分 布 对 SEP 在 行星 际 空间 中 传播 的 影响 . 通过 求解 聚焦 传输 方程 ,我 们 的 

= 模拟 结果 表明 :(1) 快 太阳 风 条 件 下 ,绝热 冷却 效应 项 发 挥 了 更 大 的 作用 ,使 粒子 能 量 衰减 的 更 快 ,而 慢 太 阳 风 对 粒 
子 的 通 量变 化 没有 显著 影响 ;(2) 加 入 观测 的 磁场 数据 时 ,粒子 的 全 向 通 量 齐 o. 

在 : 通 量 峰 值 推迟 到 达 出现 多 峰 结 构 .各 向 异性 也 发 生 一 些 改变 . 分析 表明 真实 磁场 的 极 性 对 粒子 在 行星 际 空 

> 中 传播 有 着 重要 的 影 We]. 
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Abstract Solar energetic particles (SEPs) pose one of the most serious hazards to spacecraft 
systems and constrain human activities in space. Thus, it is of importance to forecast SEP 
events. Several theories and numerical models are applied to simulate SEP events. Each model 


makes some assumptions to simplify the complex acceleration and transportation processes within 
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such events. In general, SEP will interact with ambient solar wind and background magnetic field 


during transportation. It is recognized that interplanetary transport effects must be taken into 


account at any analysis of SEP propagation. In the previous simulation, it always assumed Parker 


magnetic field and fixed solar wind speed as the input parameters. However. these assumptions 


are too simple when compared with the real conditions. In order to get better results, it is 


necessary to use more accurate background conditions. Recently. we change the fixed solar wind 


speed into spatial-dependent speed profile based on Parker's theory. and replace the Parker 


magnetic field with another Parker-like magnetic field based on in situ data at 1 AU. By solving 


the focused transport equation with simulation of time-backward stochastic processes method, our 


results show that: (1) Under fast solar wind speed assumption, it is clear that the omnidirectional 


flux decreases faster than that for the situation with slow solar wind speed in the decay phase. We 


suggest that it is due to the adiabatic cooling effect. Fast solar wind speed has a significant effect 


on the adiabatic cooling. which leads the SEPs to lose energy more quickly during transportation. 


However, slow solar wind speed has less impact on the time profiles of SEP flux and anisotropy. 


We also compare the time profiles of SEP event observed at different observatories and energies, 


the results remain the same as previous; (2) When applying in situ data of magnetic field 


observed by WIND during different Carrington Rotations, the omnidirectional flux time profiles 


vary greatly, and the main results are as followings: the peak flux appears to be delayed. multi- 


peak occur, anisotropy also has some differences. We think it results from the magnetic field 


polarity, which affects the pitch angle, and, furthermore, modulates the momentum. The 


characteristics are similar in solar minimum and solar maximum, while the peaks seem to be more 


when solar activity is active. We conclude that the real magnetic field polarity may exert a 


significant influence during the propagation of SEP. In the future, we will try to use the real-time 


background conditions which obtain from MHD models in our simulations. in order to make a 


thorough study of the SEP propagation. 
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太阳 高 能 粒子 (Solar Energetic Particles. 简称 
SEP) 事件 是 能 量 粒 子 的 通 量 突然 增强 的 事件 , 美 
国 国 家 海洋 和 大 气管 理 局 (NOAA) 将 SEP. 事件 定 
义 为 能 量 大 于 10 MeV 的 粒子 在 连续 15 min 以 上 
的 时 间 内 数目 超过 10 pfuCcm ^ e s ' e sr OR 
发 性 事件 (Rodriguez-Gasen, 2011). SEP. 事件 主要 
包括 三 种 类 型 (Kallenrode,2003): 与 太阳 光斑 爆发 
相关 联 的 脉冲 型 事件 ,与 日 晃 物 质 抛射 驱动 的 激 波 
相关 联 的 组 变型 事件 ,以 及 同时 具有 组 变型 和 脉冲 
型 事件 特征 的 混合 型 事件 . SEP. 事件 会 对 卫星 和 罕 
航 员 构成 严重 的 威胁 ,如 能 准确 预报 SEP 事件 ,人 
们 便 可 以 对 卫星 和 星 载 设备 采取 必要 的 防护 措施 ， 
同时 也 可 以 为 航天 器 的 故障 分 析 和 航天 员 的 安全 防 


lin. 
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件 产生 和 发 展 的 物理 机 制 非 常 复 杂 , 目 前 还 没有 被 
完全 解释 清楚 ,预报 SEP 事件 的 能 力也 不 尽 如 人 意 
(Lario.2005). 因此 ,正确 理解 和 恰当 地 动力 学 描述 
这 些 机 制 ,对 推进 和 提高 空间 天 气 预 报 能 力 是 非常 
重要 的 ( 魏 稳 稳 等 ,2015). 

SEP 事件 中 的 高 能 粒子 在 不 同 的 时 间 尺 度 内 
TOR DEI Tes: CME 了 驱动 的 激 波 加 速 , 并 且 在 被 
扰动 后 的 行星 际 磁场 和 太阳 风 中 传输 . 近 几 十 年 里 ， 
很 多 研究 者 探究 了 太阳 风 背 景 场 对 高 能 粒子 在 行星 
际 中 传播 的 影响 . Lario 等 (2003,2008) 基 于 Ulysses 
以 及 ACE 卫星 的 观测 ,对 太阳 活动 高 年 里 较 低能 
段 的 高 能 粒子 事件 进行 了 分 析 , 认 为 行星 际 磁 力 线 
结构 会 影响 高 能 粒子 的 经 纬度 分 布 . 并 重点 分 析 了 
大 尺度 的 行星 际 结构 对 2004 年 9 月 SEP :R TERIS 
响 , 认 为 其 强度 和 各 向 异性 时 间 剖 面 是 由 一 个 共 转 
的 低 密 度 .低速 、 低 质子 8 值 的 太阳 风流 决定 . 大 尺 


护 提 供 科学 依据 , 尽 可 能 地 降低 损失 . 然而 SEP 事 


度 太 阳 风 结构 对 高 能 粒子 传播 的 影响 决定 不 同位 置 
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的 宇宙 飞船 观测 到 的 SEP 事件 的 性 质 ,不 同 于 低 8 
值 区 域 有 利于 高 能 粒子 的 散射 ,被 压缩 的 磁场 区 域 
会 阻碍 高 能 粒子 的 自由 流动 . 磁 镜 效应 与 散射 过 程 
能 够 对 高 能 粒子 进行 有 效 的 约束 ,这 些 结构 相对 于 
观测 者 和 粒子 源 项 的 位 置 决定 SEP 事件 的 特征 
(Lario et al. ,2008). 另外 ,Shen 等 (2008) 的 研究 也 
表明 磁 云 这 种 复杂 结构 也 能 对 高 能 粒子 进行 有 效 的 
约束 并 显著 影响 高 能 粒子 的 通 量 . 

近年 来 一 些 研 究 者 采用 聚焦 传输 方程 模拟 粒子 
在 行星 际 空间 中 的 传播 过 程 , 聚 焦 传 输 方程 包括 了 
许多 重要 的 粒子 传输 机 制 : 沿 磁 力 线 的 流动 .太阳风 
对 流 .投掷 角 扩 散 、 磁 场 聚焦 效应 .垂直 扩散 、 绝 热 冷 
却 效 应 等 .模拟 结果 (Ruffolo, 1995; Qin et al., 
2006; Zhang et al. ,2009) 发现 ,在 聚焦 传播 方程 中 
引入 绝热 冷却 效应 项 ,将 增加 粒子 强度 的 衰减 率 . 
Lario 等 (Lario et al. ,1997) 研 究 了 太阳 风 对 能 量 为 
0. 1—20 MeV 的 质子 在 行星 际 空间 传播 的 影响 ,发 
现 太阳 风 对 流 能 使 粒子 更 早 地 到 达观 测 者 ,绝热 冷 
却 效应 则 导致 粒子 通 量 减 小 . 

SEP 事件 中 的 高 能 粒子 在 被 扰动 后 的 行星 际 
磁场 和 太阳 风 中 传输 ,这 一 过 程 紧 紧 依赖 于 太阳 风 
和 背景 场 (Barouch and Burlaga,1976;Lario et al. . 
1997; Malandraki et al. ,2007). 以 往 粒 子 传输 模型 
中 大 多 采用 理想 的 Parker 螺旋 磁场 和 恒定 的 太阳 
风速 度 , 背景 场 的 设置 过 于 简单 (Heras et al. , 
1995; Li et al. ,2003; Zhang et al. ,2009; Wang et 
al. ,2012). 因此 ,建立 更 加 接近 物理 真实 的 背景 场 
来 模拟 粒子 在 三 维 行星 际 空间 中 的 传播 是 提高 
SEP 事件 模拟 结果 不 可 或 缺 的 部 分 ,本 文 将 围绕 背 
景 场 模拟 展开 讨论 ,主要 研究 太阳 风速 度 以 及 真实 
磁场 分 布 对 SEP 在 行星 际 空间 中 传播 的 影响 . 


2 模型 简介 


SEP 在 日 球 层 中 的 传播 是 SEP 研究 中 的 一 个 
核心 问题 ,SEP 的 空间 分 布 等 都 与 之 密切 相关 . 研 
究 高 能 粒子 在 行星 际 空 间 中 的 传输 时 需要 求解 聚焦 
传输 方程 . 聚焦 传输 方程 的 本 质 是 Fokker-Plank 7r 
程 ,用 于 描述 粒子 分 布 函 数 f 在 相 空 间 的 演化 . 有 
些 学 者 (Kallenrode and Wibberenz, 1997; Ruffolo, 
1995; Lario et al. ,1997) 通 过 有 限 差 分 格式 对 聚焦 
传输 方程 进行 离散 求解 . Qin 等 (2006) 和 Zhang 等 
(2009) 开 发 并 发 展 了 一 种 计算 SEP 在 三 维 日 球 层 
磁场 中 传播 的 模式 ,该 模式 采用 后 向 随机 微分 方法 


解 聚 焦 传 输 方 程 .用 于 模拟 研究 观测 者 处 于 行星 际 
空间 不 同 的 经 度 、 纬 度 以 及 径 向 距离 时 ,SEP 通 量 
和 各 向 异性 随时 间 的 演化 过 程 . 

模型 采用 三 维 的 聚焦 传输 方程 来 研究 粒子 的 传 
播 过 程 ,其 形式 如 下 : 


of . 3p Əf 
TD —N. eV A oJ. 
at KI Í + (vub 十 Vas) f T Ou 


[oye 


aO — P ov sy 3b, |2 
at FE (Ve Va 7 3b: V V.) 


9n 
l—£ (y... — b. VV.O Eid bb V pE 
xi 2 toV sw $ sw ctu ; JA 


=0, (D 
其 中 f£ pe P9 ds PRI. x， 是 垂直 扩散 系数 ， 
L HBI B As UA Leo 是 粒子 的 速度 ,b 是 沿 着 磁 
场 方向 的 单位 矢量 , V。 是 太阳 风 的 速度 , D,, 为 投 
掷 角 扩散 系数 , Ls = (b * VInB) ! 为 磁场 聚焦 尺度 ， 
p 是 粒子 的 动量 . 

方程 中 的 绝热 冷却 项 的 表达 式 如 下 : 


dp — 1 一 好 E 3 
i ME (Ve Vo bb V V. +p bb; Va |, 


(2) 


RIARI y 随时 间 演 化 的 表达 式 为 
dy (Ly) | pu. — 
| 2 


Dc... —3bb.V 
dt 2Ls (V e V, —3bb:V Va). 


(3) 
我 们 根据 Qin 等 和 Zhang 4$ (Zhang, 1999; Qin 
et al. 2006; Zhang et al. ,2009) 采 用 的 后 向 随机 微 
分 方法 来 求解 聚焦 传输 方程 . 方程 (1) 可 以 改写 为 五 
个 一 维 的 后 向 随机 微分 方程 : 


dx(s) = 2k, dw(s) + (Voki— uu GO b — Vw)ds 


dCs) = J/2maxCD,, ,0) dvw(s) 4 E | (1 — pu G»v 


E BER 


y DEEE S Art 
Wa ; EGD .六 LANI 


BESTE 
di EA EM 


X pls)ds. 
求解 随机 微分 方程 采用 后 向 时 间 的 技巧 , 即 让 
虚拟 粒子 从 指定 的 地 点 ro 出 发 ,以 一 定 的 初始 动量 
户 , 从 时 刻 = to 开始 沿 后 向 时 间 运 动 ,直到 粒子 退 
出 边界 .我们 假设 所 有 的 高 能 粒子 都 在 太阳 附近 注 
入 ,因此 ,粒子 分 布 函 数 在 源 区 满足 边界 条 件 : 
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fiGegsp.D | = Pe CD aca. (4) 


其 中 , a(8,0) 是 一 个 关于 日 j 和 日 面 纬度 0 
的 函数 ,用 来 描述 SEP 源 区 强度 的 空间 分 布 ,p 是 
源 区 粒子 动量 ,T. 和 T p 
面 的 上 升 和 衰减 时 间 尺 度 . 我 们 以 下 的 模拟 事例 中 
均 假 设 源 区 里 面 的 粒子 分 布 均匀 一 致 , 即 a(y,0) — 1. 

方程 的 解 即 为 所 有 数值 模拟 过 程 在 退出 点 上 的 
平均 值 : 


enno 2 
其 中 , x。 ,ne,p。 表示 虚拟 粒子 退出 边界 时 的 参量 , 
是 粒子 退出 边界 的 时 刻 , N 是 虚拟 粒子 的 总 量 . 
粒子 全 向 通 量 (Omnidirectional flux) 的 计算 公式 


A MED ui 0 fi D 


各 向 异性 (Anisotropy) 指数 A 的 定义 式 为 A = 


fes Dm en "VN 
ee PATHS 


的 参数 见 表 1CZhang et al. ,2009) ,我 们 的 模拟 过 程 
不 考虑 垂直 扩散 对 粒子 传播 的 影响 ,将 垂直 扩散 系 
数 设 置 为 0, BẸ K1 一 0. 


R1 模型 中 参数 设 定 


Table 1 Model parameters used in the calculations 

参数 公式 初始 值 设 定 
内 边界 r—R,. R,—0.05 AU 
外 边界 r=Ro Ro=50 AU 
太阳 自转 us 2 . 

Q = 2x/T T=27. 27 
频率 ar K 
太阳 风 Vsw = Vswer Vsw = 400 km * s^! 
B= Bo Zoa BoQsing, 

磁场 pe UY. Bo=3.54 nT * AU? 
投掷 角 D,, = Do(1— gl)vp7h Do 为 常数 


b=2/3, h—0.2, q—5/3 


扩散 系数 a 
(a) = (os? = (B/ B.) 


X GT |a| Tk) 


3. 模拟 结果 及 分 析 


3.1 太阳 风速 度 对 粒子 传播 过 程 的 影响 

根据 物理 性 质 和 起 源 区 域 的 不 同 , 太阳风 一 般 
分 为 快 太阳 风 (v 三 550 km * s) RUE EH (Cos 
450 km * s ! (Wang et al. ,2000). 我 们 分 别 计算 
了 快 . 慢 太阳 风 对 粒子 通 量 的 影响 , 快 太阳 风 对 不 同 
能 量 质 子 的 传播 过 程 以 及 对 不 同 观测 位 置 的 通 量 剖 


面 的 影响 . 并 且 统 计 分 析 了 快 、 慢 太阳 风 条 件 下 统计 
到 的 粒子 的 动量 分 布 特征 . 模型 中 磁场 采用 Parker 


螺旋 磁场 . 

为 了 研究 发 展 的 太阳 风速 度 对 粒子 在 行星 际 空 
间 中 传播 的 影响 ,我 们 将 模型 中 沿 径 向 恒定 不 变 的 
太阳 风速 度 V = Vawe (Vo = 400 km。s 2 变 为 


Parker 太阳 风速 度 解 ( 由 (6) 式 解 出 ),Parker 太阳 
风速 度 解 的 具体 数值 实现 方法 详 见 文献 Parker 
(1958). 从 (6) 式 中 可 以 看 出 ,Parker 太阳 风速 度 解 
FRERE MS 日 心 距离 + 相关. 


w = Vswer, 
Hp V, mM 
1 Vi o4 YDo ( Vo a | GMe 
2 (y — 1) V r r 
Ta YDo GM, 
= 一 Vo 4 s 6 
22 l (y— Dp To n 


1 给 出 了 太阳 风速 度 接近 慢 太 阳 风 的 模拟 结 
TR. 图 la 给 出 了 该 事例 中 太阳 风速 度 随 日 心 距离 x 
的 变化 曲线 ,太阳 风速 度 在 1AU 处 达到 383 km * s". 
以 下 图 形 中 标注 为 慢 太 阳 风 的 都 是 基于 图 1a 中 的 
Parker 慢 太 阳 风 速度 解 . 图 1b 给 出 了 不 同 太阳 风 
速度 条 件 下 100 MeV 能 量 的 质子 在 1AU 赤道 处 的 
全 向 通 量 对 比 图 . 黑色 实 线 代表 固定 速度 400 km * s ' 
的 模拟 结果 ,红色 虚线 代表 固定 速度 383 km s ^! 
的 模拟 结果 ,绿色 虚线 代表 Parker 太阳 风速 度 解 接 
近 慢 太阳 风速 度 下 的 模拟 结果 . 从 图 1b 中 可 以 看 出 
这 三 条 通 量 曲线 非常 接近 ,说 明 慢 太 阳 风 相对 固定 
的 太阳 风速 度 而 言 ,对 粒子 的 通 量变 化 影响 不 大 . 

2 给 出 了 太阳 风速 度 接近 快 太 阳 风 的 模拟 结 
果 . 图 2a 中 Parker 太阳 风速 度 解 在 1AU 处 达到 
644 km * s .以 下 图 形 中 标注 为 快 太阳 风 的 都 是 
基于 图 2a 中 的 Parker 快 太阳 风速 度 解 .图 2b 给 出 
了 不 同 太 阳 风 速度 条 件 下 100 MeV 能 量 的 质子 在 
1AU 赤道 处 的 全 向 通 量 对 比 图 . 黑色 实 线 代 表 固 定 
速度 400 km。s 的 模拟 结果 ,红色 虚线 代表 固定 
速度 644 km* s 的 模拟 结果 ,绿色 虚线 代表 Parker 
太阳 风速 度 解 接近 快 太阳 风速 度 下 的 模拟 结果 . 从 
2b 中 可 以 看 出 ,固定 太阳 风速 度 变 大 时 ,将 增加 
粒子 通 量 曲 线 的 衰减 率 ,而 Parker 快 太阳 风 条 件 
下 ,粒子 通 量 曲 线 的 衰减 率 进一步 增加 .通过 分 析 我 
们 认为 这 是 由 于 快 太阳 风 条 件 下 ,绝热 冷却 效应 项 
发 挥 了 更 大 的 作用 . 在 Parker 太阳 风速 度 解 的 模型 
中 ,绝热 冷却 项 dp/dt 的 表达 式 如 下 : 
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图 1 太阳 风速 度 接近 慢 太 阳 风 的 模拟 结果 
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Ca). 太阳 风速 度 随 日 心 距离 的 变化 ; (b) 不 同 太阳 风 条 件 下 能 量 为 100 MeV 的 质子 的 全 向 通 量 对 比 . 


Fig.1 Simulation results under slow solar wind speed assumption 


(a) Slow solar wind speed profile; (b) Comparison of 100 MeV protons flux between different solar wind speed. 
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图 2 太阳 风速 度 接近 快 太阳 风 的 模拟 结果 
度 随 日 心 距离 的 变化 ; (b) 不 同 太阳 风 条 件 下 能 量 为 100 MeV 的 质子 的 全 向 通 量 对 比 . 


(a) 太 阳 风 速 
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Fig.2 Simulation results under fast solar wind speed assumption 


(a) Fast solar wind speed profile; (b) Comparison of 100 MeV protons flux between different solar wind speed. 


2 = -a[o , Vm — bb: V.) but bbs Vo | p 


- p sin! g+ p cosp Sr T | 


Vs dea 
r| 2 


第 一 项 ， Le. Vu —bbIVVS)— Le [(2 — sin!) V. /r ^- sin'BdV,, /dr.]; 


第 二 项 H 7s bb :VV、， 


从 式 (7) 可 以 看 出 ,粒子 的 动量 随时 间 的 函数 是 
衰减 的 ,绝热 冷却 效应 会 使 得 观测 到 的 粒子 比 它 在 
源 区 释放 时 的 能 量 要 小 . 相 比 于 沿 径 向 恒定 不 变 的 
太阳 风速 度 ,在 Parker 太阳 风速 度 解 的 条 件 下 , 绝 
热 冷却 项 dp/dt 右 侧 表达 式 中 的 方 括号 里 多 了 


sing m Ru cos? g T T 项 ,这 两 项 在 快 


= p’ sin? BV /rp Tor 


| 
慢 太 阳 风 条 件 下 均 为 正 值 . 在 快 太阳 风 条 件 下 ,不 仅 
方 插 号 里 的 数值 增 大 ,V 也 超过 400 km * s ,使 
得 表达 式 右 侧 显著 变 小 ,导致 绝热 冷却 效应 相对 速 

度 固 定 为 400 km * s 的 太阳 风 明 显 增强 . 而 在 慢 
太阳 风 条 件 下 ,虽然 方 括号 里 的 数值 增 大 ,但 是 Vs。 

小 于 400 km* s “ ,使 得 表达 式 右 侧 变化 不 大 ,导致 
绝热 冷却 效应 相对 而 言 没 有 显著 改变 . 男 外 ,SEP 
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的 谱 指数 是 小 于 零 的 ,所 以 对 于 某 一 特定 能 段 的 粒 
子 而 言 , 随 着 时 间 推 移 ,观测 到 的 粒子 通 量 会 越 来 越 
少 . 因此 ,在 快 太阳 风 条 件 下 ,粒子 的 通 量 如 图 2 所 
示 训 减 的 更 快 ,粒子 全 向 通 量 在 后 期 可 达到 一 个 量 
级 的 差异 . 而 慢 太 阳 风 条 件 下 则 如 图 1 所 示 没 有 显 
著 变化 . 
由 于 在 SEP 的 初始 相 时 期 粒子 密度 的 梯度 大 ， 
扩散 效应 比 绝热 冷却 效应 更 重要 ,所 以 在 SEP. 的 初 
始 相 可 以 忽略 绝热 效应 ,在 SEP 的 衰减 时 期 ,粒子 
密度 空间 梯度 小 ,绝热 冷却 效应 将 起 重要 的 作用 ( 涂 
传 论 等 ,1988). 在 Parker 太阳 风速 度 解 的 条 件 下 ， 
对 于 某 一 固定 的 径 向 距离 7, 太阳 风速 度 增 大 时 . 4H 
应 的 Vw/r 和 dV/dr 也 随 之 增 大 ,第 一 项 和 第 二 
项 与 太阳 风速 度 的 大 小 呈 线 性 关系 ,因此 ,太阳 风速 
度 越 高 ,会 造成 SEP 向 外 对 流 的 越 多 ,绝热 能 量 损 
失 也 随 之 增 大 ,绝热 效应 也 就 越 显 著 . 

我 们 又 对 快慢 太阳 风速 度 条 件 下 统计 到 的 粒 
子 的 动量 进行 了 统计 分 析 . 图 3 给 出 了 观测 到 的 粒 
子 的 动量 分 布 图 . 黑色 实 线 ,红色 虚线 分 别 代表 
Parker 太阳 风速 度 解 接近 快 太阳 风速 度 . 慢 太阳 风 


表 2 不 同 太 阳 风 条 件 下 的 统计 结果 
Table 2 Statistical results under different solar 


wind speed conditions 


Din Pmax Ap ”统计 到 的 粒子 总 数 
V fast 0. 448 1. 507 0. 106 201463 
V sow 0. 447 1. 008 0. 056 215300 


(1) 从 所 统计 到 粒子 的 最 低 动量 来 看 ,它们 都 
比 目 标 粒 子 对 应 的 动量 高 ,说 明 粒 子 在 运动 过 程 中 
都 有 能 量 衰减 ; 

(2) 从 所 统计 到 粒子 的 最 高 动量 来 看 , 快 太阳 
风 条 件 下 粒子 能 量 的 最 大 值 更 大 ,说 明快 太阳 风 条 
件 下 粒子 能 量 的 衰减 可 以 更 加 显著 ; 

(3) 投放 相同 的 粒子 数 (投放 1X10" 个 粒子 ) 
时 , 快 太阳 风 条 件 下 统计 到 的 粒子 数 小 于 慢 太 阳 风 
条 件 下 的 粒子 数 ,也 验证 了 上 述 结论 ; 

(4) 从 粒子 分 布 来 看 , 越 靠 近 目 标 粒 子 能 量 段 
的 粒子 数 越 多 , 且 呈 现 出 单调 递减 的 趋势 ,说 明 大 部 
分 统计 到 的 粒子 是 能 量 衰减 较 少 的 粒子 . 

此 外 ,我 们 又 模拟 了 快 太阳 风 对 不 同 能 量 的 粒 


速度 时 的 统计 结果 . 我 们 将 观测 到 的 粒子 动量 的 最 
大 值 和 最 小 值 (p,。。& pw ) 等 分 为 10 个 区 间 Ap 
0.1(pmax 一 pmin) ,然后 统计 每 个 动量 区 间 内 的 粒子 
数 占 总 的 粒子 数 的 百分比 . 表 2 为 不 同 太 阳 风 条 件 
下 的 统计 结果 ,我 们 设置 的 目标 粒子 的 动量 为 


子 和 不 同 观 测 位 置 处 的 通 量 剖 面 的 影响 . 图 4a、4b 
分 别 给 出 了 不 同 太 阳 风 速度 条 件 下 100 MeV 能 量 
的 质子 在 观测 者 位 于 0. 5AU 赤道 处 以 及 2AU 、 纬 
BE 80 处 的 全 向 通 量 对 比 .图 4c 给 出 了 不 同 太阳 风 
速度 条 件 下 10 MeV 能 量 的 质子 在 LAU 赤道 处 的 全 向 


5 —0. 445. 其 中 目标 粒子 的 动量 定义 为 观测 者 处 最 
终 观测 到 的 粒子 动量 ,通过 能 量 计算 可 以 得 到 ,该 能 
是 人 为 设 定 的 ,是 为 了 便于 统计 分 析 某 一 个 具体 
量 段 内 的 粒子 的 传播 特征 .分布 特征 等 . 

从 图 3 和 表 2 的 结果 可 以 看 出 : 
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图 3 快慢 太阳 风 条 件 下 统计 到 的 粒子 动量 分 布 


Fig. 3 Momentum distributions of observed particles 


通 量 对 比 图 . 黑色 实 线 代 表 固 定 速度 为 400 km * s f] 
模拟 结果 ,红色 虚线 代表 Parker 太阳 风速 度 解 接近 
快 太阳 风速 度 下 的 模拟 结果 . 从 结果 中 可 以 得 出 , 快 
太阳 风 对 不 同 能 量 段 的 粒子 在 不 同 观测 位 置 处 的 通 
量 剖 面 有 着 相似 的 影响 . 
3.2 用 实时 的 磁场 观测 数据 

以 往 粒 子 传输 模型 中 大 多 采用 Parker 螺旋 人 磁 
场 ,Parker 螺旋 磁场 是 一 个 稳 态 模型 ,在 太阳 风 是 
定 态 球 对 称 的 假定 下 ,根据 麦克 斯 韦 方程 组 、 磁 冻结 
效应 ,以 及 在 源 表面 CR. ,90.,$) 处 径 向 的 边界 条 件 
可 以 推导 出 Parker 螺旋 磁场 模型 ,具体 形式 为 (Parker， 


1958) 

B= B, (Ey e p, (Ey Bur EUM, 
其 中 0 为 太阳 自转 频率 ,V, 为 常数 太阳 风速 度 . 对 
于 距离 太阳 很 远 的 地 方 7 传 R.. 根据 磁 通 量 守恒 定 
律 ,模型 中 用 到 的 是 简化 的 Parker 磁场 ,如 表 1 中 


所 示 : B = Be, ono, 


r 


V $? 其 中 BSQ—3.54 nT 是 
r SW 
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图 4 快 太阳 风 条 件 下 观测 者 处 于 不 同位 置 时 能 量 为 100 MeV 的 质子 的 全 向 通 量 对 比 ((a) 和 (b)) 
以 及 能 量 为 10 MeV 的 质子 的 全 向 通 量 对 比 (c) 


Fig.4 Comparisons of 100 MeV protons flux at different locations (panels (a) and (b)) and energies (panel (c)) 
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图 5 (a) 不 同行 星际 磁场 下 ,对 比 100 MeV 质子 的 通 量 ; Co 各 向 异性 ; (b) 根据 WIND 飞船 数据 得 到 的 


CR2056-CR2057 内 磁场 B, 分 量 图 . 太阳 坐标 为 (X,Y) = (rcosg,rsing) = (0,0)， 黑 色 圆 圈 为 地 球 


Fig.5 Distribution of 100 MeV proton omnidirectional flux (a) and anisotropy (c) at different IMF conditions 


during CR 2056-2057. (b) is distribution of B, based on Wind spacecraft data. The Sun is at (X,Y) = (rcosy, 
rsing) = (0,0) and the Earth is at the solid black circle 


由 LAU 处 的 磁场 佑 算得 到 . 间 中 的 传输 也 可 能 发 挥 重要 作用 .因此 ,为 了 构造 更 
然而 ,这 种 简化 了 的 Parker 磁场 忽略 了 磁场 大 加 真实 的 背景 磁场 ,我 们 使 用 了 WIND 卫星 实时 的 
小 和 极 性 的 变化 ,但 是 这 些 因 素 对 粒子 在 行星 际 空 “ ”磁场 观测 数据 来 进行 背景 磁场 的 设置 . Florens 等 
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《2007) 推 导出 了 可 以 根据 卫星 实时 观测 数据 求解 磁 
场 三 分 量 (B, .B,.D,) 的 类 Parker 磁场 解 , 如 下 式 
(8) 所 示 ,我 们 根据 该 方法 将 WIND 卫星 的 观测 数 
据 拓展 到 全 空间 . 

QG — R, )sing, 


B= Bo (È) e -Bgo (È 1j v 5 
(8) 
假设 测试 粒子 从 太阳 源 表 面 (R.,y.) 离开 ,到 
达 空 间 任 意 位 置 Cr. 40 . 并 且 太 阳 风 在 每 一 个 卡 林 
顿 周期 内 是 不 变 的 , 则 可 以 根据 宇宙 飞船 的 观测 数 
据 估 算 R, 处 太阳 风 的 参量 . 参数 n R4 ,上 办， 观测 
时 间 z, KELo. to 3- T ]BJ JT A NER] 由 以 下 公式 
计算 得 到 : 
4. = QG —4) —QG, 
$ = 34,0 Q(r — RO/V.. 


模型 中 使 用 的 是 WIND 飞船 一 分 钟 精度 的 磁场 数 
据 ,我 们 分 别 计算 了 处 于 太阳 活动 低 年 以 及 太阳 活 
动 高 年 的 不 同 卡 林 顿 自转 周期 (简称 CR) 内 SEP 的 
全 向 通 量 和 各 向 异性 的 时 间 剂 面 . 
太阳 活动 低 年 的 模拟 结果 如 图 5 所 示 . 其 中 ,图 
5b 是 根据 WIND 飞船 在 2007 年 5 月 3 日 至 2007 
年 5 月 29 日 处 于 CR2056 和 CR2057 内 的 观测 数据 
得 到 的 黄道 面 内 的 磁场 图 . 为 了 更 好 的 可 视 化 效果 ， 
pie 日 心 距离 7 处 的 磁场 进行 了 标准 化 处 
理 , 具 体 标 准 化 处 理 方法 为 :将 同一 条 Parker 螺 线 
上 ,任意 日 心 距离 x 处 的 磁场 , 均 赋 为 该 条 螺 线 在 
LAU 处 的 磁场 数值 ,然后 进行 可 视 化 作 图 . 太阳 坐 
标 为 (X,Y) = (rcosg,rsing) = (0,0) , 地 球 处 在 黑 
色 圆 圈 处 ,并 按照 顺 时 针 方向 进行 背景 磁场 的 设置 . 
两 条 黑 线 的 夹 角 为 统计 到 的 粒子 所 在 源 区 的 大 致 分 
布 范围 . 图 5a 中 红色 虚线 为 在 此 CR 期 间 内 根据 类 
Parker 磁场 模拟 得 到 的 粒子 全 向 通 量 剖 面 ,黑色 实 
线 为 原 Parker 磁场 下 的 模拟 结果 ,图 5c 为 相应 的 
各 癌 异 性 剖面 . 相 比 之 下 ,我 们 发 现在 类 Parker 磁 
场 下 模拟 得 到 的 粒子 通 量 峰值 会 到 达 的 晚 一 些 , 全 
向 通 量 剖 面 会 出 现 多 个 峰值 并 旦 在 衰减 相 下 降 的 更 
快 , 而 且 一 阶 各 向 异性 也 有 明显 的 不 同 . 
由 式 子 (10) 一 (12) 可 知 , 当 加 入 磁场 的 观测 数 
据 后 , cosp 和 sing 的 大 小 不 会 改变 ,但 是 它们 的 符 
号 会 随 着 磁场 的 极 性 而 发 生 相 应 的 改变 ,从 而 会 影 
WERA u 值 的 变化 ,并 最 终 对 动量 p 的 值 产生 影 
响 . 因此 ,我 们 认为 粒子 通 量 峰 值 推迟 到 达 , 多 峰 结 
构 以 及 各 向 异性 的 差异 主要 是 由 于 加 入 的 真实 磁场 


R/V w> 
(9) 


具有 不 同 的 极 性 所 致 . 

针对 通 量 峰值 推迟 到 达 , 我 们 比较 了 Parker 螺 
旋 磁 场 和 类 Parker 螺旋 磁场 下 最 大 峰值 附近 的 粒 
子 初始 能 量 的 大 小 ,发现 Parker 磁场 下 粒子 初始 能 
量 的 平均 值 大 约 为 0. 463, 28 Parker 磁场 下 粒子 初 
台 能 量 的 平均 值 大 约 为 0. 628, 而 我 们 统计 的 目标 
粒子 能 量 为 0. 445. 由 于 能 量 越 高 ,衰减 到 目标 能 量 
所 需要 的 时 间 越 长 ,所 以 在 类 Parker 磁场 中 由 于 磁 
场 极 性 对 粒子 能 量 的 调制 ,它们 到 达观 测 者 所 需 的 
时 间 会 更 久 , 从 而 导致 峰值 推迟 到 达 . 

针对 出 现 的 多 峰 结构 ,我 们 在 图 6 中 给 出 了 第 
二 个 较 大 峰值 到 达 前 后 不 同时 间 段 内 统计 到 的 粒子 
初始 动量 p 的 分 布 范 围 . 黑色 实 线 、 红 色 虚 线 、 绿 色 
虚线 分 别 代 表 第 二 个 峰值 到 达 前 两 天 ,第 二 个 峰值 
期 间 以 及 第 二 个 峰值 过 后 两 天 内 的 粒子 动量 分 布 . 
从 图 6 可 以 看 出 ， RR edic Mida 
显 差异 .在 第 二 个 峰 出 现 的 时 间 段 内 (红色 虚线 )， 
hs i xb ON p? RS 
函数 ) 的 计算 ,我 们 发 现 这 个 峰值 对 该 期 间 的 粒子 通 
量 的 贡献 起 决定 性 作用 ,并 且 使 得 该 通 量 超 过 了 前 
后 两 天 的 通 量 大 小 . 通常 情况 下 随 着 时 间 的 推移 ,在 
衰减 相 不 同时 间 段 内 统计 到 的 大 部 分 粒子 的 动量 会 
不 断 增 大 , 即 类 似 图 6 中 黑色 实 线 的 峰值 会 不 断后 
移 , 从 而 使 得 全 向 通 量 不 断 减 小 . 然而 ,同样 ,由 于 磁 
场 极 性 对 粒子 能 量 的 调制 ,使 得 粒子 动量 出 现 了 额 
外 的 峰值 并 发 挥 了 主要 作用 ,从 而 导致 出 现 了 多 峰 
结构 . 
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图 6  CR2056-CR2057 模拟 得 到 的 不 同时 间 段 内 
到 达观 测 者 的 粒子 初始 动量 p 的 分 布 范 围 
Fig. 6 Initial momentum distributions of observed particles 


in different intervals during CR2056-CR2057 
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不 同 于 Parker 磁场 下 的 粒子 各 向 异性 ,类 Parker 
磁场 下 的 各 向 异性 在 上 升 相位 时 小 于 零 ,而 在 衰减 相 
期 间 , 各 向 异性 在 模拟 误差 范围 内 会 在 零 附近 上 下 浮 


动 . 由 于 观测 点 的 磁场 极 性 保持 不 变 , 则 根据 各 向 异性 
Do te ie 
人 


以 及 y= 二 p*b/p 可 知 , 当 A 为 正 ( 负 ) 时 ,统计 到 的 
主要 是 朝 磁 场 方向 同 ( 反 ) 向 运动 的 粒子 . 因此 ,在 上 
升 相位 期 间 , 朝 磁场 相反 方向 运动 Ca — — D. 的 粒子 
占 主 导 . 而 在 衰减 相 期 间 , 由 于 磁场 极 性 的 影响 以 及 
部 分 粒子 投掷 角 方 向 的 改变 ,会 使 得 各 向 异性 在 误 
差 范围 内 出 现 波动 . 

太阳 活动 高 年 的 模拟 结果 ,我 们 以 CR 1969 为 
例 展示 在 了 图 7 中 . 类 似 图 5 的 格式 ,图 Ta 给 出 了 
模拟 得 到 的 粒子 全 向 通 量 剖 面 ,图 7c 为 各 向 异性 剖 
面 ,图 7b 给 出 了 WIND 飞船 在 该 周期 内 所 观测 到 
的 磁场 . 其 中 两 条 黑色 线 同样 代表 统计 到 的 粒子 所 
在 源 区 的 大 致 范围 . 类 似 太阳 活动 小 年 的 情况 ,从 图 
7a 和 图 7c 中 可 以 看 出 ,粒子 通 量 峰 值 也 比 原 Parker 
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峰 结 构 ,而 各 向 异性 在 误差 范围 内 并 没有 太 大 的 差异 . 
同样 地 ,我 们 计算 得 到 类 Parker 磁场 下 最 大 峰 
值 附近 所 统计 到 的 粒子 初始 能 量 的 平均 值 大 约 为 
0.562, 要 大 于 Parker 磁场 下 粒子 能 量 的 平均 值 0. 463， 
因此 粒子 在 类 Parker 螺 线 磁场 里 运动 到 达观 测 者 
所 需 的 时 间 更 久 ,导致 峰值 推迟 到 达 . 但 是 , 相 比 于 
太阳 活动 小 年 粒子 的 平均 初始 能 量 为 0. 628 ,该 周期 
内 统计 到 的 粒子 初始 能 量 要 更 小 一 些 , 所 以 相 比 之 
下 粒子 的 全 向 通 量 剖面 峰值 到 达 的 也 更 早 一 些 . 至 
于 更 小 的 粒子 初始 能 量 , 可 能 与 太阳 活动 高 年 更 复 
杂 的 磁场 极 性 变化 对 粒子 能 量 的 调制 有 关 . 
此 外 ,与 太阳 活动 小 年 相 比 ,在 太阳 活动 高 年 里 
模拟 得 到 的 粒子 的 全 向 通 量 剖 面 会 出 现 更 多 的 峰 
值 . 图 8 给 出 了 第 二 ,三 个 峰值 到 达 时 以 及 它们 到 达 
之 前 一 段 内 观测 到 的 粒子 初始 动量 p 的 分 布 范围 . 
从 图 8 可 以 看 出 ,在 这 两 个 峰值 到 达 之 前 ,如 黑 线 所 
示 ,粒子 初始 动量 大 约 以 0. 8 为 峰值 ,而 且 粒 子 的 动 
量 分 布 范围 更 加 广泛 .在 红线 所 示 的 第 二 个 峰值 到 
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Ca) 不 同行 星际 磁场 下 ,对 比 100 MeV 质子 的 通 量 ; (c) 各 向 异性 ; 
(bo 根据 WIND 飞船 数据 得 到 的 CR1969 的 磁场 B. 分 量 

Fig.7 Distribution of 100 MeV proton omnidirectional flux (a) and anisotropy (c) at different IMF conditions 

during CR 1969. (b) is distribution of B, based on Wind spacecraft data 


图 7 
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图 8 CR1969 模拟 得 到 的 不 同时 间 段 内 到 达观 测 者 的 
粒子 初始 动量 p 的 分 布 范 围 


Fig.8 Initial momentum distributions of observed 


particles in different intervals during CR1969 


达 期 间 , 粒 子 的 初始 动量 更 加 集中 地 分 布 在 以 0. 6 左 
右 为 峰值 的 一 个 小 范围 内 ,因此 根据 通 量 计算 公式 
可 知 该 段 时 间 内 粒子 的 全 向 通 量 更 大 . 而 在 绿 线 所 
示 的 第 三 个 峰值 到 达 期 间 ,粒子 的 初始 动量 仍 以 0. 6 
左右 为 峰值 ,但 是 它们 占 该 时 间 段 内 所 统计 到 总 粒 
子 数 的 比例 更 小 ,而 且 来 源 更 加 广泛 ,这 就 是 为 什么 
第 三 个 峰值 要 比 第 二 个 峰值 低 的 原因 . 类 似 于 太阳 
活动 小 年 的 情况 ,在 太阳 活动 高 年 里 磁场 对 粒子 能 
量 的 调制 会 导致 多 峰 结构 ,而 这 种 更 多 峰值 结构 的 
出 现 , 是 由 于 磁场 的 极 性 出 现 了 更 多 的 变化 ,对 粒子 
能 量 的 调制 更 加 复杂 所 致 ，. 

从 图 7c 中 可 以 看 出 ,类 Parker 磁场 下 的 各 向 


和 背景 磁场 对 粒子 在 行星 际 空 
拟 结果 表明 : 

COD 将 固定 的 太阳 风速 度 改 为 变化 的 Parker 
太阳 风速 度 解 ,从 模拟 结果 可 以 看 出 慢 太 阳 风 对 粒 
子 的 通 量变 化 没有 显著 影响 ,而 当 太 阳 风 速度 加 快 
时 ,粒子 通 量 剖 面 在 后 期 下 降 的 更 快 . 这 是 因为 快 太 
阳 风 条 件 下 ,绝热 冷却 效应 更 加 显著 ,从 而 使 粒子 的 
能 量 衰减 的 更 快 . 

(2) 相对 于 原 Parker 磁场 模型 ,加 入 观测 的 磁 
场 数据 时 粒子 的 全 问 通 量 剖 面 发 生 了 比较 明显 的 变 
化 : 通 量 峰 值 推迟 到 达 、 出 现 多 峰 结构 ,各 向 异性 也 
发 生 一 些 改变 . 通过 分 析 ,我 们 认为 加 入 观测 的 磁场 
数据 时 ,磁场 有 了 极 性 ,粒子 的 投掷 角 会 随 着 磁场 的 
极 性 而 发 生 相 应 改变 ,进而 也 会 对 粒子 的 能 量 进 行 
调制 ,从 而 导致 模拟 结果 出 现 了 上 述 的 变化 . 无论 在 
太阳 活动 低 年 还 是 高 年 ,粒子 全 向 通 量 齐 面 和 各 向 
异性 剖面 都 伴随 着 这 些 改变 ,而 太阳 活动 高 年 更 加 
复杂 的 变化 可 能 是 由 于 磁场 的 极 性 更 加 复杂 所 和 致 . 

我 们 的 模拟 结果 说 明太 阳 风 速度 以 及 背景 磁场 
对 粒子 在 行星 际 空 间 中 的 传播 有 着 重要 的 影响 . 以 
往 粒子 传输 模型 中 大 多 采用 恒定 的 太阳 风速 度 和 
Parker 螺旋 磁场 ,背景 场 的 设置 太 过 理想 ,这 可 能 
导致 无 法 准确 模拟 太阳 高 能 粒子 在 行星 际 空间 中 的 
演化 过 程 . 因此 ,设置 更 加 真实 的 背景 场 对 提高 我 们 
的 模拟 结果 有 着 重要 的 意义 . 未 来 的 研究 工作 中 我 
们 将 尝试 加 入 由 MHD 模拟 (Feng et al. , 2010,2012; 
Shen et al. , 2011,2014) 得 到 的 实时 变化 的 太阳 风速 
度 和 三 维 磁场 分 量 , 建 立 更 加 完整 .更 加 接近 物理 真 
实 的 太阳 风 背 景 场 ,以 期 得 到 更 真实 的 粒子 在 三 维 


间 中 传播 的 影响 , 模 


异性 在 上 升 相 位 时 大 于 零 ,说 明 朝 磁场 相同 方向 运 
动 ( — 1) 的 粒子 占 主导 . 同样 地 ,在 衰减 相 期 间 ， 
由 于 磁场 极 性 的 影响 以 及 粒子 投掷 角 方 向 的 改变 ， 
会 使 得 各 向 异性 在 误差 范围 内 出 现 一 些 波动 . 

从 模拟 结果 可 以 看 出 ,太阳 活动 高 年 和 低 年 的 
结果 虽然 都 包含 通 量 峰 值 推迟 到 达 、 出 现 多 峰 结 构 、 
各 向 异性 发 生 改 变 等 特点 ,但 也 存在 差异 . 由 于 太阳 
活动 高 年 磁场 比较 复杂 ,导致 通 量 剖 面 出 现 更 多 的 
峰 , 各 向 异性 剖面 也 存在 差异 .但 二 者 的 模拟 结果 都 
充分 说 明了 真实 磁场 极 性 对 粒子 在 行星 际 空间 中 的 
传播 有 着 重要 的 影响 . 


4 结论 


本 文通 过 求解 聚焦 传输 方程 研究 了 太阳 风速 度 


行星 际 空间 的 传播 物理 过 程 , 同 时 可 以 研究 日 冕 物 
质 抛 射 相互 作用 等 对 高 能 粒子 事件 形成 和 传播 的 影 
i] Gopalswamy et al. , 2002; Shen et al. ,2008; Li 
et al. , 2012; Shen et al. ,2013). 
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